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Properties of jets measured from tracks in proton-proton collisions at center-of-mass 

energy ^/s = 7 TeV with the ATLAS detector 
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Jets are identified and their properties studied in center-of-mass energy ^/s = 7 TeV proton- 
proton collisions at the Large Hadron Collider using charged particles measured by the ATLAS 
inner detector. Events are selected using a minimum bias trigger, allowing jets at very low trans- 
verse momentum to be observed and their characteristics in the transition to high-momentum fully 
perturbative jets to be studied. Jets are reconstructed using the anti-fct algorithm applied to charged 
particles with two radius parameter choices, 0.4 and 0.6. An inclusive charged jet transverse mo- 
mentum cross section measurement from 4 GeV to 100 GeV is shown for four ranges in rapidity 
extending to 1.9 and corrected to charged particle-level truth jets. The transverse momenta and 
longitudinal momentum fractions of charged particles within jets are measured, along with the 
charged particle multiplicity and the particle density as a function of radial distance from the jet 
axis. Comparison of the data with the theoretical models implemented in existing tunings of Monte 
Carlo event generators indicates reasonable overall agreement between data and Monte Carlo. These 
comparisons are sensitive to Monte Carlo parton showering, hadronization, and soft physics models. 

PACS numbers: 13.85.Hd,13.87.Ce,13.87.Fh 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) [1, 2] provides an 
excellent description of the kinematic distribution of high 
transverse momentum jets in proton-proton collisions, 
but does not give straightforward predictions for the 
properties of particles within these jets or for the prop- 
erties of low-momentum jets. These quantities may be 
predicted by Monte Carlo event generators, whose re- 
sults are dependent on phenomenological models of par- 
ton showering, hadronization, and soft (i.e. low momen- 
tum transfer) physics. These models have free parame- 
ters that must in turn be tuned to data. 

In this work the properties of low-momentum jets are 
measured from charged particle tracks in the ATLAS 
2010 data at center-of-mass energy \/s = 7 TeV. These 
properties will be compared to a range of Monte Carlo 
tunes derived from previously measured data, allowing 
the study of the transition between the MC generators' 
separately-implemented models of soft strong interac- 
tions and perturbative QCD. This measurement is com- 
plementary to others that are used for comparison with 
models of soft QCD, such as inclusive charged particle 
kinematic distributions [3] and underlying event proper- 
ties [4] , and may shed light on the limited ability of exist- 
ing models to describe both types of measurements simul- 
taneously with a single set of tunable parameters [5] . At 
higher momenta, track-based jet measurements provide 
a complement to calorimeter-based measurements [6, 7], 
with results independent of calorimeter calibrations, se- 
lections, and uncertainties. 

Jets have previously been measured using charged par- 
ticles in proton-proton collisions at the CERN Intersect- 
ing Storage Rings [8] , and in proton-antiproton collisions 



at the CERN Super Proton Synchrotron [9] and Fermi- 
lab Tevatron [10]. The momentum fraction of charged 
particles with respect to calorimeter jets has also been 
measured [11, 12]. Related fragmentation measurements 
in electron-positron and electron-proton collisions are re- 
viewed in Ref. [13]. 

In this work, jets are reconstructed using the anti-/ct 
algorithm [14], using a radius parameter R of 0.4 or 0.6, 
applied to measured charged particles with transverse 
momentum px > 300 MeV; jet four-momenta are deter- 
mined by adding constituent four-vectors. The distribu- 
tion of jet momenta, jet track multiplicity, and kinematic 
properties of tracks within these track-based jets are cor- 
rected back to truth-level charged particle jets, which are 
defined to be the jets obtained when the same algorithm 
is applied to all primary charged particles emerging from 
the proton-proton collision with the same px cut. 

Five quantities are measured for charged particle jets: 
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The inclusive charged particle jet cross section ajet is 
measured in bins oipTjet for four rapidity (y) [15] ranges: 
0-0.5, 0.5-1.0, 1.0-1.5, 1.5-1.9, where the last range is 
narrower than the others to ensure that measured jets 
arc fully contained in the tracking acceptance of the 
ATLAS detector. The distribution of charged particle 
multiplicities per jet N'j'j^^ is measured for each of these 
ranges, separated further into five jet transverse mo- 
mentum ranges: 4-6 GeV, 6-10 GeV, 10-15 GeV, 15- 
24 GeV, 24-40 GeV. The longitudinal momentum frac- 
tion z of charged particles in these jets is measured in 
the same rapidity and transverse momentum ranges. The 



variable z, also known as the fragmentation variable, is 
defined for each particle in a jet by: 



Pch ■ Pjet 
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where Pch is the momentum of the charged particle and 
Pjet is the momentum of the jet that contains it. The z 
distribution presented here differs from the usual defini- 
tion (in [11], for example), which would include neutral 
particles and low-momentum charged particles in the to- 
tal jet momentum. The variable p^'' is the momentum 
of charged particles in a jet transverse to that jet's axis: 
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Finally, the density of charged particles in (j>-y space, 
Pch{r), is measured as a function of the radial distance r 
of charged particles from the axis of the jet that contains 
them, where r is given by 



ch - (pjetY + [Uch - VjetY 



SO that 
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Note that this is a particle number density, rather than 
the related energy density variable used for calorimeter- 
based jet shape measurements. The variables Pch{r) and 
pJ^^ are measured in the same rapidity and jet momentum 
ranges as z. 



II. THE ATLAS DETECTOR 

The ATLAS detector [16] at the Large Hadron Colhder 
(LHC) [17] covers almost the whole solid angle around the 
collision point with layers of tracking detectors, calorime- 
ters and muon chambers. It has been designed to study 
a wide range of physics topics at LHC energies. For the 
measurements presented in this paper, the tracking de- 
vices and the trigger system are used. 

The ATLAS inner detector (ID) has full coverage in 
the pseudorapidity range \r]\ < 2.5. It consists of a sil- 
icon pixel detector (Pixel), a silicon microstrip detector 
(SCT) and a transition radiation tracker (TRT), the last 
of which only covers \ri\ < 2.0. These detectors cover 
a sensitive radial distance from the interaction point of 
50-150 mm, 299-560 mm and 563-1066 mm, respectively, 
and are immersed in a 2T axial magnetic field. The 
inner-detector barrel (end-cap) parts consist of 3 (2x3) 
Pixel layers, 4 (2x9) double- layers of single-sided silicon 
microstrips with a 40 mrad stereo angle, and 73 (2x 160) 



layers of TRT straws. Typical position resolutions are 
10, 17 and 130 fim. for the R-cp coordinate of the pixel de- 
tector, SCT, and TRT respectively; the resolution of the 
second measured coordinate is 115 (580) ^m for the pixel 
detector (SCT). A track from a charged particle travers- 
ing the barrel detector would typically have 11 silicon 
hits (3 pixel clusters and 8 strip clusters) and more than 
30 straw hits. 

The ATLAS detector has a three-level trigger system: 
Level 1 (LI), Level 2 (L2) and Event Filter (EF). The LI 
triggers use custom fast electronics; the L2 and EF fur- 
ther refine selections in software, with L2 using a subset 
of the event information for rapid processing and EF re- 
constructing complete events. For this measurement, the 
trigger relies on the LI signals from the Beam Pickup 
Timing devices (BPTX) and the Minimum Bias Trigger 
ScintiUators (MBTS). The BPTX stations are composed 
of electrostatic button pick-up detectors attached to the 
beam pipe at ±175 m from the center of the ATLAS de- 
tector. The coincidence of the BPTX signal between 
the two sides of the detector is used to determine when 
bunches are colliding in the center of the detector. The 
MBTS arc mounted on each side of the detector at Z == 
±3.56 m. They are segmented into eight sectors in az- 
imuth and two rings in pseudorapidity (2.09 < |?7| < 2.82 
and 2.82 < \r]\ < 3.84). Data were collected for this anal- 
ysis using a trigger requiring a BPTX coincidence and 
MBTS trigger signals. The MBTS trigger used for this 
paper is configured to require one hit above threshold 
from either side of the detector, referred to as a single- 
arm trigger. 



III. MONTE CARLO SAMPLES 

The corrections used in this analysis, from detector- 
level track jets to truth-level charged particle jets, are 
derived from a Monte Carlo (MC) simulated sample us- 
ing the Pythia 6.421 event generator program [18] with 
the ATLAS AMBTl tune, whose parameters are cho- 
sen based on single track distributions in ATLAS mini- 
mum bias data [3]. In order to derive corrections, these 
events are then passed through the ATLAS detector sim- 
ulation [19], based on GEANT 4 [20]. Large simulated 
samples of the ATLAS MC09 [21] and Perugia 2010 [22] 
tunes of Pythia 6 are also used. 

Additional Pythia 6.421 tunings are used at the 
generator-level for comparison to data, and in order to 
explore the sensitivity of the measurement to underlying 
event, fragmentation, and hadronization parameters, as 
described in Section VII. The Perugia family of tunes [22] 
uses the CTEQ5L [23] parton distributions; there is a 
central value ( "Perugia 0" ) and several variants that at- 
tempt to bracket the possible changes in tune parameters 
that are permitted in fitting existing data. Two vari- 
ants, "Perugia HARD" and "Perugia SOFT", are cre- 
ated by changing the initial-state radiation cutoff scale 
so that the pcrturbative contribution is different, while 



other parameters are adjusted to bring the tuned distri- 
butions back into agreement. Perugia HARD (SOFT) 
has more (fewer) jets per event with higher (lower) av- 
erage momentum, and a higher (lower) average charged 
multiplicity. The Perugia 2010 variant is an adjustment 
of Perugia that improves the description of jet shapes 
and adjusts hadronization parameters for better consis- 
tency with LEP data, given that the original values for 
these parameters were based on Q^ shower ordering; in 
particular, it has an adjusted yield of strange particles. 
Perugia NOCR ("no color reconnection" ) represents an 
effort to reproduce the same data with no color recon- 
nection model used; it is used in computing systematic 
uncertainties, but because it could not be made to fit all 
the original tuning data it is not included on comparison 
plots. All Pythia tunes described thus far use px-ordered 
showering; the DW tune is also used in order to include 
the impact of Q^ ordering. 

Predictions from separately-implemented generators 
are used for additional studies and comparisons. Phojet 
1.12.1.3 [24], which relies on Pythia 6.115 for the frag- 
mentation of partons, is used with its default settings. A 
Herwig-I— I- 2.4.2 [25] sample is used with its default set- 
tings, along with an additional sample tuned on 7 TeV 
underlying event data with Herwig-|~+ 2.5.1 (UE7) [26]. 
A Pythia 8.145 [27] sample is used with the 4C tune [28]. 

Single- and double-diffractive events are included in 
the Phojet and Pythia 8 samples, with cross sections rel- 
ative to the non-diffractive events as indicated by those 
generators. HerwigH — V does not include diffractive mod- 
eling, while Pythia 6 diffractive modeling produces a neg- 
ligible yield of jets as defined in this analysis. 



IV. DATA SELECTION 

This measurement uses a sample of early ATLAS data 
at center-of-mass energy ^/s = 7 TeV for which the mini- 
mum bias trigger was minimally prescaled, corresponding 
to a total integrated luminosity of 799 /xb~^ and a peak 
luminosity of 6.6 x 10^^ cm~^s~^. Events from colliding 
proton bunches are selected if the MBTS recorded one 
or more counters above threshold on cither side. They 
are further required to have a primary vertex [29] re- 
constructed using beam-spot information [30]. Events 
with additional reconstructed primary vertices are re- 
jected. The average number of collisions per bunch cross- 
ing /x depends on luminosity; the highest-luminosity data 
collected have ^ ^ 0.14, but over half the data have 
A* < 0.01. 



A. Track reconstruction 

Tracks are reconstructed using the ATLAS primary 
silicon-seeded tracking algorithm, with a configuration 
similar to that described in the 900 GeV and 7 TeV min- 
imum bias measurements [3, 31]. In order to select good- 



quality tracks emerging from the primary vertex while 
maintaining a high efficiency, each track is then required 
to have: 

• transverse momentum px > 300 MeV; 

• pseudorapidity \ri\ < 2.5; 

• transverse impact parameter with respect to the 
primary vertex \do\ < 1.5mm (0.2mm) for tracks 
with pt less than (greater than) 10 GeV; 

• longitudinal impact parameter with respect to the 
primary vertex zq satisfying jzosin^j < 1.5mm; 

• if a signal (or hit) is expected in the innermost pixel 
detector layer (i.e. if the extrapolated track passes 
through a section of that layer with functioning in- 
strumentation) , then such a hit is required, with 
one pixel hit in any layer required otherwise; 

• at least 6 SCT hits. 



B. Selection efficiencies 

The efficiency for triggering given the presence of one 
or more jets is determined from data using a random 
subset of all events with a BPTX coincidence. In such 
events for which an offline track jet (as defined in Sec- 
tion V) is present, the fraction of events for which the 
MBTS single-arm trigger is also passed is determined. It 
is found to be negligibly different from 100% for all events 
containing jets as defined in this analysis. 

The efficiency for primary vertex-finding given the 
presence of one or more jets is determined from data 
by removing the track selection cuts that use impact pa- 
rameter with respect to the primary vertex described in 
Section IV A. Other cuts are kept as described, and a re- 
quirement is added that the transverse impact parameter 
with respect to the beam spot satisfies |do| < 4.0mm. 
Jets are then reconstructed from tracks satisfying this 
new selection, and the fraction of events containing jets 
that also have a primary vertex reconstructed is deter- 
mined. It is found to be negligibly different from 100% 
for all events containing jets as defined in this analysis. 

The inefficiency introduced by the requirement that 
no additional primary vertices be present arises from two 
causes: 

• Events may be lost due to multiple collisions in a 
single bunch-crossing. The rate at which this oc- 
curs is calculated as a function of /z [32] , and cor- 
responds to a 3.3% correction to the cross section 
over the entire dataset. 

• Events may be lost because multiple vertices are 
incorrectly reconstructed in bunch crossings with 
only a single collision. The rate at which this oc- 
curs is determined from simulated events, and pa- 
rameterized as a function of the number of selected 
tracks in the event. It corresponds to a cross sec- 
tion correction of roughly 2%, depending weakly on 
jet momentum. 



The efficiencies for track reconstruction and selection 
are derived from simulated events, whose tracking prop- 
erties have been shown to be in excellent agreement with 
the data [3]. The uncertainties on these efficiencies, and 
their impact on the present results, are discussed in Sec- 
tion VII. 



V. JET DEFINITION 

Reconstructed track jets are produced by applying the 
anti-fct algorithm with radius parameter i? = 0.4 or 0.6 to 
the selected tracks; the pion mass is assumed for all tracks 
in the application of the algorithm. After the algorithm 
is applied, track jets with px > 4 GeV and |y| < 1.9 
are accepted, including those with only one constituent 
track. The rapidity cut is chosen to ensure that jets are 
fully contained within the tracking acceptance \ri\ < 2.5 
in the R ~ 0.6 case; for the narrower R = 0.4 jets, the 
same rapidity cut is used for consistency. 

Truth-level charged particle jets similarly have the 
anti-fc( algorithm applied to primary charged simulated 
particles with px > 300 MeV, using each particle's 
true mass in the application of the algorithm. Primary 
charged particles are defined as charged particles with a 
mean lifetime r > 0.3 x 10~^'^s, which are produced in 
the primary collision or from subsequent decays of par- 
ticles with a shorter lifetime. Thus the charged decay 
products of Kg particles arc not included. Charged par- 
ticle jets are required to meet the pt and y requirements 
given for track jets above. 



VI. CORRECTION PROCEDURE 

After all distributions are measured for reconstruction- 
level track jets, each distribution is corrected to truth- 
level charged particle jets. The corrections are derived 
from the AMBTl sample discussed in Section III, and 
account for tracking efficiency and momentum resolution. 

Reconstructed jets arc binned simultaneously in px, 
constituent multiphcity (Nj^^), and rapidity (y), so that 
a three-dimensional distribution is produced for the pur- 
pose of applying corrections. Similarly, for each per-track 
variable z, p^\ and r, each track in each jet is binned 
in the variable, the parent jet momentum, and parent jet 
rapidity. For both the jet-level quantities and each track- 
level variable, corrections are applied simultaneously in 
the three binned variables using the Baycsian Iterative 
Unfolding algorithm [33]. 

The corrections for this algorithm are based on a re- 
sponse matrix derived from MC events, which encapsu- 
lates the probability for a charged particle jet with a par- 
ticular {pT,Nj^^,y) to be reconstructed in each possible 
px, Nj^^, and y bin. A reconstructed track jet is de- 
fined to be matched to a truth-level charged particle jet 
if it is within AR < 0.2 (0.3) for jet radius 0.4 (0.6). 
These values are chosen to avoid matching ambiguity 



while still giving good efficiency, which rises from roughly 
55% at the very lowest momenta to greater than 95% for 
PTjet > 10 GcV. In the case of the z response matrix, 
binned in {z,pTjet,yjet), if a truth jet is unmatched then 
all its tracks are counted as lost to inefficiency. If it is 
matched, then each of its constituent charged particles 
are matched to the closest track within AR < 0.04; if 
there is no such track, then that particular particle is 
counted as lost to inefficiency. The pJj?' and r matrices 
are filled in the same manner as the z. 

The Bayesian Iterative Unfolding algorithm takes both 
inefficiencies and resolution effects into account, as 
"missed" entries and bin-to-bin transfers respectively; 
the probability of these occurrences is determined from 
the central AMBTl MC sample. Before the Bayesian It- 
erative Unfolding algorithm can be applied, a correction 
must be made for unmatched reconstructed jets, and for 
unmatched reconstructed tracks in the z, p'^\ and r dis- 
tributions. This correction is determined from the central 
MC sample by the ratio between the reconstructed dis- 
tributions and the matched reconstructed distributions. 

Three iterations of the algorithm are used. This is val- 
idated by applying the unfolding to ensembles of data- 
sized MC samples, and observing the number of itera- 
tions required to obtain convergence. Statistical uncer- 
tainties on the unfolded variables are determined by a toy 
histogram- variation method, which produces 100 pseudo- 
experiment distributions by varying the contents of each 
input bin randomly, then calculating the covariancc ma- 
trix from the variation in the output of the unfolding. 
This error calculation is also validated using ensembles 
of data-sized MC samples (with full detector simulation) 
by considering the distribution of pulls: 
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where Ncorr is the corrected number of events in a given 
bin, Ntruth is the true number of events in that bin, and 
Gcorr is the Statistical error for that bin as reported by the 
unfolding procedure. The RMS of the pull distribution 
should be consistent with one. In some bins it is not, 
so the statistical errors are scaled appropriately. This 
scaling is at most roughly 1.5 except for a few outlying 
bins in the multiplicity and pT^^ distributions that have 
very low statistics. 

The ensemble studies also reveal a bias in the unfolding 
algorithm calculation in closure tests, which is roughly 
2% in the jet px distributions. The possibility of an ad- 
ditional bias from changing the prior distribution used as 
input to the unfolding algorithm was investigated. After 
three iterations, this bias vanishes. 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are computed by determin- 
ing the impact of a given effect on reconstruction-level 



distributions, then propagating it through the unfolding 
procedure to determine the final impact on the measure- 
ment. A summary of all systematic uncertainties for se- 
lected bins is shown in Table I for the cross section and 
p?P^' distributions, and in Tabic II for the multiplicity, z, 
and p{r) distributions. The uncertainties given in each 
column of the tables are: 

Tracking efficiency and resolution: The systematic 
uncertainties associated with tracking efficiency 
and resolution are computed by applying efficiency 
scaling and resolution smearing to single tracks in 
generated events. The efficiencies and resolutions 
are derived from fully-simulated events and param- 
eterized by a track's 1/pt, Vj ^^^ the AR to the 
nearest other track. These efficiency uncertain- 
ties are determined by adjusting the derived single- 
track efficiencies in accordance with uncertainties 
determined from data [3] , for which uncertainties in 
the material description of the ATLAS inner detec- 
tor dominate. The resulting distributions are com- 
pared with those found for the baseline resolution, 
with the increase (decrease) taken as the upper 
(lower) systematic uncertainty. Similarly, the un- 
certainties in tracking resolution arc investigated by 
applying an additional Gaussian resolution smear- 
ing derived from studies of high- momentum muons; 
the resulting change is taken as a symmetric uncer- 
tainty. 

Monte Carlo: The efficiency scaling and resolution 
smearing are similarly used for uncertainties asso- 
ciated with the particular event generator used to 
compute the unfolding corrections. For each tune 
described in Section III, an ensemble of 50 samples 
with the same number of events as the data has the 
smearing applied, and the average unfolded distri- 
bution is compared to one produced similarly (with 
196 data-sized samples) for the central AMBTl MC 
sample. The largest increase (decrease) is taken 
as the upper (lower) systematic uncertainty. Al- 
though the response matrix primarily models the 
impact of detector efficiency and resolution, truth- 
level details can impact the corrections in several 
ways. Different jet or track momentum distribu- 
tions change the population within bins; differences 
in the amount or distribution of activity in the un- 
derlying event can significantly change the track 
momentum and radius distribution, especially in 
the low-momentum bins, because these tracks are 
not correlated with the jet direction and so may ap- 
pear at large radii. Differing strangeness fractions 
can change the distribution of long-lifetime tracks 
that decay and produce kinks in tracks, leading to 
momentum mismeasurements and/or loss of tracks 
due to failed hit requirements. The variation over 
tunes also accounts for PDF uncertainties, because 
MC properties are tuned with a particular PDF 
and several different PDF's are used. 



High-pT Tracks: Despite the tightening of the do cut 
for high-momentum tracks, 0.6% (8.3%) of selected 
tracks in the simulation with momentum above 
10 (40) GeV do not have a matching truth par- 
ticle with momentum within 50%. These high- 
momentum mismeasurcd tracks arc primarily asso- 
ciated with wide-angle scatterings in the material of 
the inner detector that create the appearance of a 
single, straight track [3]. Were the fraction of such 
tracks similar in data, their impact on the measure- 
ment would be accounted for in the unfolding pro- 
cedure. However, the data have a larger fraction of 
high-pT tracks failing the do cut (2.3%) than does 
the MC (1.6%). A systematic uncertainty on the 
number of mismeasured tracks in data, as a func- 
tion of track px, is computed and propagated to 
assess its impact on jets. The upper bound of this 
uncertainty assumes that all such extra tracks are 
well-measured, but did not pass the dp cut due to 
improperly modeled resolution. The lower bound 
is based on the assumption that the increase in re- 
jected tracks corresponds to a proportional increase 
in accepted mismeasured tracks. These uncertain- 
ties are then propagated to each measured jet bin, 
with the scale factor determined by the fraction 
of jets in a given jet bin with their leading track 
in each momentum range. Thus the correction is 
largest for high jet momentum and low number of 
particles per jet, because these jets have the highest 
leading track momenta. This uncertainty on the re- 
constructed data is then used to scale the measured 
distributions, with the unfolding applied in each 
case and compared to the unfolded central value. 
The resulting differences in each measured bin give 
the final uncertainty due to high- momentum mis- 
measured tracks. This correction is also applied to 
the z, p^'' and r distributions. Each bin is scaled in 
proportion to the fraction of tracks in that bin that 
come from a jet whose leading track is in a given 
range, in proportion to the uncertainty on tracks 
in that range. The corrections arc only significant 
for high-z bins in high-pT jets. The corrections are 
then propagated through the unfolding as in the 
jet-based distributions, and the result compared to 
the central value. 

Unmatched jets and tracks: In order to assign an 
uncertainty to the correction for unmatched re- 
constructed jets and tracks, the unfolding is re- 
peated with a correction determined from a fully- 
simulated Perugia 2010 in place of the baseline AT- 
LAS AMBTl. The difference in the unfolding out- 
put between this and the baseline configuration is 
taken to be the uncertainty on the correction. The 
uncertainty is symmetrized by taking the maximum 
deviation. 

Split vertex: The data are corrected for event rejection 
due to misreconstruction of extra primary vertices 



TABLE I. 
An overall 



Summary of systematic uncertainties for selected bins in selected cross section and p^' distributions, for R = 0.6. 
luminosity uncertainty of 3.4% is not shown, and applies to the cross section only. 



Distribution 


Bin 

[GeV] 


Track 
eff. 


Track 
res. 


Monte 
Carlo 


High-pT 
Tracks 


Unmatched 
Jets/Tracks 


Split 
Vertex 


Clo- 
sure 


Total 


d'^Tjct 


4-5 


+3.2% 
-3.3% 


±0.07% 


+0.07% 
-1.2% 


< 0.005% 


±0.21% 


±1.8% 


±2.0% 


+4.2% 
-4.4% 


dp-Tj^tdyjct 

U/jetKO.5 


20-22 


+6.6% 
-6.3% 


±0.34% 


+2.2% 
-0.69% 


+0.10% 
-0.14% 


±0.05% 


±2.0% 


±0.54% 


+7.3% 
-6.7% 


40-45 


+6.8% 
-6.9% 


±0.39% 


+ 1.4% 
-2.3% 


+ 1.1% 
-1.4% 


±0.01% 


±2.1% 


±0.50% 


+7.3% 
-7.7% 


80-90 


+7.1% 
-6.9% 


±0.34% 


+7.1% 
-7.5% 


+5.5% 
-9.1% 


< 0.005% 


±2.3% 


±0.78% 


+ 12% 
-14% 


d^CTjet 


4-5 


+5.4% 
-5.3% 


±0.02% 


+3.0% 
-2.1% 


< 0.005% 


±0.19% 


±1.8% 


±5.2% 


+8.3% 
-7.9% 


1.5<|!;jct|<1.9 


20-22 


+ 10% 
-9.5% 


±0.02% 


+2.6% 
-3.2% 


+0.10% 
-0.15% 


±0.09% 


±2.0% 


±0.96% 


+ 11% 
-10% 


40-45 


+ 11% 
-10% 


±0.30% 


+2.9% 
-1.5% 


+ 1.1% 
-1.3% 


±0.07% 


±2.1% 


±0.95% 


±11% 


80-90 


+ 12% 
-11% 


±1.3% 


+7.3% 
-12% 


+7.3% 
-8.2% 


±9.3% 


±2.0% 


±6.9% 


+20% 
-21% 


1 dN,^ 

lH3ctl<l-9 
4GeV<pTjet<6GcV 


0.0-0.05 


+0.44% 
-0.32% 


±0.18% 


+ 16% 
-4.8% 


< 0.005% 


±2.7% 


< 0.005% 


±0.17% 


+ 16% 
-5.6% 


0.5-0.55 


+0.30% 
-0.20% 


< 0.005% 


+3.0% 
-1.7% 


< 0.005% 


±3.2% 


±0.01% 


±0.09% 


+4.4% 
-3.6% 


0.85-0.9 


+ 1.0% 
-0.03% 


±1.2% 


+ 10% 
-6.5% 


< 0.005% 


±2.9% 


±0.02% 


±2.1% 


+ 11% 
-7.5% 


1 diV.fe 
iVjet dp-' 

iHjetKl-S 

24 GGV<pTj^t <40 GoV 


0.0-0.05 


+ 1.2% 
-1.1% 


±0.10% 


+3.4% 
-3.0% 


+0.21% 
-0.17% 


±0.58% 


±0.05% 


±0.05% 


+3.6% 
-3.2% 


0.5-0.55 


±1.2% 


±0.01% 


+3.3% 
-2.4% 


+0.38% 
-0.29% 


±0.40% 


±0.03% 


±0.01% 


+3.5% 
-2.7% 


0.85-0.9 


+0.95% 
-1.0% 


±0.06% 


+3.3% 
-2.6% 


+0.32% 
-0.25% 


±0.57% 


±0.05% 


±0.09% 


+3.5% 
-2.9% 


3.0-3.5 


+2.6% 
-2.4% 


±0.56% 


+8.7% 
-1.1% 


+0.05% 
-0.04% 


±0.27% 


±0.14% 


±0.48% 


+9.1% 
-2.7% 



as a function of the number of seleeted tracks in 
the event, as discussed in Section IV B. A weight 
is applied on an event- by- event basis based on the 
probability of that event being rejected. As the cor- 
rection is derived entirely from simulation, the full 
value of the correction is taken as a (symmetrized) 
uncertainty. 

Closure: As discussed in Section VI, closure tests reveal 
a bias in the correction procedure. This is taken as 
a bin-by-bin systematic uncertainty, which is sym- 
metrized. 

The total uncertainty on the cross section is dom- 
inated by tracking efficiency, with mismeasured high- 
momentum tracks also playing a role for the highest- 
momentum jets. Tracking efficiency uncertainties play 
a similar role for the multiplicity distributions, especially 
for low-multiplicity jets, with Monte Carlo uncertainties 
also making large contributions. For other distributions, 
Monte Carlo uncertainties are dominant, with closure un- 
certainties making large contributions in extremal bins; 
tracking efficiency uncertainties mostly cancel because 
these distributions are normalized by the number of jets. 



VIII. RESULTS AND DISCUSSION 

A selection of the distributions measured in this anal- 
ysis appears in Figures 1-6. Other rapidity ranges may 



be found in Ref. [34]. They are compared to the MC 
distributions described in Section III. 



A. Charged particle jet cross section 

Cross sections as a function of jet px are shown in Fig- 
ure 1 . The simulated cross sections shown for comparison 
are scaled to the data, using the scale factor S defined 
by: 



Q __ total /total 
*-* — ^data I ^MC 



(7) 



where 



total 



100 GoV 



4 GoV 



dV 



jet 



dpTjet / dyjet-, , 

/_1.9 dpTjetdyjet 



(8) 



The scale factors for the various MC's are given in Ta- 
ble III. 

In all cases except IIerwig±+ 2.4.2, the scale factors S 
for R = 0.6 are larger than those for R = 0.4. Since the 
larger radius parameter results in the inclusion of more 
particles not directly associated with perturbative scat- 
tering, this implies that the models underestimate the 
contribution of the underlying event required to repro- 
duce the data. Phojet has the best agreement between 
the scale factors at i? = 0.4 and R = 0.6. The Pythia 



TABLE II. 
0.6. 



Summary of systematic uncertainties for selected bins in selected multiplicity, 2, and pchir) distributions, for R = 



Distribution 


Bin 


Track 
eff. 


Track 
res. 


Monte 
Carlo 


High-pT 
Tracks 


Unmatched 
Jets/Tracks 


Split 
Vertex 


Clo- 
sure 


Total 


1 dJVj.t 
4 GoV <pT3et<6GcV 


1 


+6.8% 
-6.3% 


±0.01% 


+7.9% 
-10% 


< 0.005% 


±0.15% 


±0.21% 


±0.25% 


+10% 
-12% 


3 


+2.8% 
-2.6% 


±0.14% 


+3.3% 

-0.27% 


< 0.005% 


±0.16% 


±0.08% 


±0.26% 


+4.3% 
-2.6% 


5 


+1.2% 
-1.3% 


±0.03% 


+0.23% 
-2.4% 


< 0.005% 


±0.03% 


±0.04% 


±0.09% 


+1.2% 
-2.7% 


9 


+ 11% 
-10% 


±0.36% 


+ 11% 
-8.4% 


< 0.005% 


±0.01% 


±0.24% 


±1.4% 


+ 15% 
-13% 


1 dN,^, 

|j/jetl<l-9 

24 GoV<pT3et <40 GcV 


1 


+ 19% 
-11% 


±3.2% 


+50% 
-8.5% 


+3.7% 
-5.5% 


±0.13% 


±0.35% 


±35% 


+64% 
-38% 


3 


+10% 
-9.9% 


±0.05% 


+3.1% 
-3.3% 


+1.7% 
-1.9% 


±0.14% 


±0.28% 


±0.77% 


±11% 


5 


+6.8% 
-6.5% 


±0.03% 


+3.3% 
-1.4% 


+0.64% 
-0.75% 


±0.06% 


±0.21% 


±0.26% 


+7.6% 
-6.7% 


9 


+1.6% 
-1.7% 


±0.14% 


+2.9% 
-0.21% 


±0.07% 


±0.02% 


±0.08% 


< 0.005% 


+3.3% 
-1.7% 


1 dN,i, 
Nj^t dz 

4GeV<pT3et<6GcV 


0.1 -0.125 


±1.7% 


±0.11% 


+ 1.5% 
-8.9% 


< 0.005% 


±2.0% 


±0.05% 


±0.03% 


+3.0% 
-9.3% 


0.5 - 0.525 


+1.2% 
-1.3% 


±0.07% 


+3.2% 
-0.54% 


< 0.005% 


±2.0% 


±0.03% 


±0.18% 


+3.9% 
-2.4% 


0.85-0.9 


+4.5% 
-4.7% 


±0.20% 


+3.7% 
-4.0% 


< 0.005% 


±0.25% 


±0.12% 


±0.33% 


+5.9% 
-6.1% 


1.0 


+7.0% 
-6.5% 


±0.08% 


+4.8% 
-16% 


< 0.005% 


±0.03% 


±0.23% 


±0.28% 


+8.5% 
-17% 


1 dN,^ 
24 GeV<PTjct <40 GcV 


0.1 -0.125 


+ 1.2% 
-1.3% 


±0.14% 


+3.0% 
-5.8% 


+0.42% 
-0.32% 


±0.01% 


±0.03% 


±0.09% 


+3.3% 
-5.9% 


0.5 - 0.525 


+3.9% 
-3.6% 


±0.29% 


+2.7% 
-0.86% 


+0.40% 
-0.65% 


±0.64% 


±0.07% 


±0.24% 


+4.8% 
-3.9% 


0.85-0.9 


+9.0% 
-9.3% 


±1.7% 


+ 1.8% 
-6.6% 


+3.3% 
-3.7% 


±0.21% 


±0.20% 


±1.3% 


+ 10% 
-12% 


0.95- 1.0 


±13% 


±1.3% 


+6.8% 
-3.5% 


+3.8% 
-5.5% 


±0.81% 


±0.30% 


±3.1% 


+ 16% 
-15% 


Pch(r) 

lH3ett<l-9 

4GoV<pTjct<6GcV 


0.0-0.01 


+7.1% 
-6.8% 


±0.12% 


+20% 
-46% 


< 0.005% 


±0.30% 


±0.22% 


±0.36% 


+21% 
-46% 


0.09-0.1 


+0.00% 
-0.06% 


±0.22% 


+2.6% 
-10% 


< 0.005% 


±2.4% 


< 0.005% 


±0.12% 


+3.5% 
-10% 


0.28-0.3 


+0.90% 
-1.1% 


±0.11% 


+2.5% 
-10% 


< 0.005% 


±2.4% 


±0.03% 


±0.03% 


+3.6% 
-11% 


0.55-0.6 


+2.4% 
-2.8% 


±0.06% 


+0.94% 
-3.1% 


< 0.005% 


±3.8% 


±0.07% 


±5.4% 


+7.1% 
-7.8% 


Pch{r) 

lH3=tl<1.9 
24 GoV<pTjct <40 GcV 


0.0-0.01 


+3.3% 
-3.4% 


±0.08% 


+4.2% 
-5.1% 


+0.48% 
-0.59% 


±0.20% 


±0.09% 


±0.20% 


+5.4% 
-6.2% 


0.09 - 0.1 


+0.56% 
-0.60% 


±0.04% 


+2.9% 
-4.6% 


+0.33% 
-0.26% 


±0.21% 


< 0.005% 


±0.05% 


+3.0% 
-4.7% 


0.28-0.3 


+2.0% 
-2.1% 


±0.07% 


+3.4% 
-8.0% 


+0.45% 
-0.35% 


±0.86% 


±0.10% 


±0.04% 


+4.1% 
-8.4% 


0.55-0.6 


+2.5% 
-2.3% 


±0.02% 


+3.8% 
-6.7% 


+0.47% 
-0.36% 


±1.3% 


±0.15% 


±0.64% 


+4.8% 
-7.2% 



Perugia tunes also agree well for the two radii, and are 
most consistent with one. 

The jet cross section distributions (Fig. 1) fall by 
6 orders of magnitude between jet momenta of 4 and 
100 GeV. The MC models considered agree broadly with 
this trend, but do not agree well in detail. By construc- 
tion of the normalization factor S, all distributions agree 
with the data in the lowest momentum bins; most also 
give qualitative agreement for the shape at the lowest 
momenta. The MC distributions diverge from the data 
in the 10-20 GeV range, with some having a harder and 
others a softer momentum dependence,. At higher px, 
many of the models' momentum dependence agrees well 
with the data. If one identifies the higher-momentum 
region as dominated by perturbative modeling and the 
low-momentum region as dominated by soft physics, this 
indicates that perturbative modeling of charged particle 
jets is in fair agreement for most of the tunes. It is the 
transition from soft physics to the perturbative region 



that is not successfully modeled. 

The Pythia models give a harder shape for the mo- 
mentum spectrum than the data below a jet px of about 
20 GeV, after which they exhibit roughly the same mo- 
mentum dependence or become slightly softer. By con- 
trast, the Phojet and Herwig++ models produce spectra 
that are softer than the data in the 10-20 GeV range 
but have relatively good shape agreement outside, al- 
though the Herwig+± 2.5.1 UE7 tune has the additional 
feature of producing too hard a spectrum at momenta 
below 10 GeV. 

The DW tune has a cross section that falls much too 
slowly in the transition region, and falls much too rapidly 
at high momentum; this effect is especially pronounced 
for R = 0.6. This suggests that the Q^-ordered showering 
used by the DW tune is less successful in modeling the 
jet momentum spectrum. 

Cross sections as a function of rapidity are shown in 
Figure 2. The MC distributions are normalized to the 
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FIG. 1. The cross section for anti-fct charged particle jets as a function of pT, with \y\ < 0.5 and radius parameter R as 
indicated. The shaded area is the total uncertainty for the corrected data distribution, excluding the overall 3.4% luminosity 
uncertainty. The data are compared to a range of theoretical results from Monte Carlo event generators, which are normalized 
to the data over the full momentum and rapidity range measured, using the scale factor S as defined in the text. The bottom 
inserts show the fractional difference between these distributions and the data. The distributions for the Perugia HARD 
(Perugia SOFT) tunc, not shown, agree qualitatively with the Perugia2010 (Perugia 0) tune. 



data in each momentum bin separately rather than to the 
scale factor S. The rapidity dependence of the cross sec- 
tion shows generally good agreement between data and 
MC. The cross section decreases only slightly with in- 
creasing rapidity at low momenta, but by a somewhat 
higher amount at higher momenta. 



B. Charged particle kinematics and multiplicity in 

jets 

The multiplicity of charged tracks per jet, for several 
momentum ranges over the full range \yjet\ < 1-9, is 
shown in Figure 3. The charged particle z, p^ , and 
Pchir) distributions are shown for same central rapidity 
range in Figures 4, 5, and 6 respectively. 

No tune describes well all of the kinematic distribu- 
tions and multiplicities of charged particles within jets. 
For z, the ATLAS AMBTl and Perugia SOFT give good 
descriptions of the data. For pch{r), AMBTl gives a good 



description. For p^'' and multiplicity, no tune correctly 
describes all data. For all distributions, Herwig-|— I- 2.4.2 
shows strong disagreement with the data, characterized 
by an excess of low-momentum particles, which is espe- 
cially pronounced for the larger jet-finding parameter R 
and at large particle r (as defined in Equation 4). Her- 
wigH — h 2.5.1 UE7 represents a significant improvement, 
so Herwig-|--|- 2.4.2 will not be discussed further. 

All remaining models give good agreement for the av- 
erage charged particle multiplicity per jet. The AMBTl 
and Perugia SOFT tunes agree well with the multiplicity 
distributions (Figure 3) for the vast majority of jets, and 
Perugia 0, Perugia 2010, and Phojet give fair agreement 
for R = 0.4, although the high-multiplicity tail in data is 
greatly underestimated by all models. Herwig-|--t- 2.5.1 
UE7 and Perugia HARD have a significant excess of low- 
multipHcity jets, while Pythia 8.145 4C and Pythia DW 
exhibit a deficit. 

The AMBTl and Perugia SOFT tunes give good agree- 
ment with the measured longitudinal momentum fraction 
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FIG. 2. The cross section for anti-fct charged particle jets as a function of rapidity, for selected momentum bins and radius 
parameter R as indicated. The shaded area is the total uncertainty for the corrected data distribution, excluding the overall 
3.4% luminosity uncertainty. The data are compared to a range of theoretical results from Monte Carlo event generators, which 
are normalized to the data separately for each momentum range. The distributions for all Perugia samples agree qualitatively, 
so only Perugia is shown. The two Herwig-f-|- tunes agree, so only Herwig-|— I- 2.4.2 is shown. 



z (Figure 4); Perugia-0 also agrees well for R = 0.4. The 
other MC's (except Herwig++ 2.4.2) agree within 30% at 
low z, but diverge more significantly at high z. Perugia 
HARD has the most significant excess of high-z particles, 
with excesses also present for Phojet and Perugia 0. The 
excess is particularly large at lower jet momenta and R 
= 0.6, suggesting that the soft physics model is charac- 
terized by fewer particles with higher momentum. By 
contrast, Pythia 8.145 4C, Pythia DW, and Herwig++ 
2.5.1 UE7 have too few high-z particles, with variations 
again larger for i? = 0.6. Pythia DW exhibits an excess at 
mid-z at low jet momenta, which is seen at progressively 
lower z values as the jet momentum increases, implying 
an excess of particles with a momentum of roughly 2 GeV 
that is not associated with jet structure. 

The transverse momentum p^'' is in fair agreement 
(~20%) at low-to-moderate values for all MC genera- 
tors except Herwig-|--|- 2.4.2. At the lowest jet momenta 
and highest measurable p^^' ' Ptiojet and Herwig-|— I- 2.5.1 
UE7 have an excess of particles, while Pythia DW has 
a deficit. At higher jet momenta, the data have more 
high-pJJ?' particles than any tune, with Perugia 2010 and 



Perugia HARD giving the closest description and Peru- 
gia SOFT the furthest. Perugia 2010 and Perugia HARD 
agree better than do the other Perugia tunes. 

The AMBTl, Pythia 8.145 4C, and Herwig-h-H 2.5.1 
UE7 tunes provide a good description of the charged par- 
ticle number density Pchir) (Figure 6) at all radii. Phojet 
and the Perugia tunes (especially SOFT) have an excess 
of particles very close to the jet axis, which is most pro- 
nounced at high jet momentum and for R ~ 0.6; Peru- 
gia 2010 agrees better in this region than do the other 
Perugia tunes. At high r, Phojet and all Pythia tunes 
except AMBTl and Perugia SOFT have too few parti- 
cles. However, the disagreement is less pronounced than 
is seen at high pT^'' , implying that high-radius soft parti- 
cles from the underlying event are better-described than 
high-radius hard radiation. 



IX. CONCLUSIONS 

A measurement is presented of the charged particle 
jet cross section as a function of transverse momentum 
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TABLE III. Scale factors S for Monte Carlo cross section nor- 
malization as defined in Equation 7. The total Monte Carlo 
cross sections are normalized to the total for the data, over 
the full momentum and rapidity ranges investigated in this 
analysis, 4 GeV < pTjet < 100 CeV and \yjet\ < 1.9. The sta- 
tistical uncertainties on these ratios are less than 0.1%. The 
systematic uncertainties are ^4gS (±5.0%) for R = 0.4 (0.6); 
these uncertainties are entirely correlated within columns, and 
largely correlated between columns. 



Tunc 


S iR = 0.4) 


S {R^ 0.6) 


AMBTl 


0.838 


0.896 


Perugia-0 


0.981 


1.087 


Perugia HARD 


0.936 


1.058 


Perugia SOFT 


0.968 


1.036 


Perugia 2010 


0.976 


1.044 


DW 


0.894 


1.045 


Herwig-I-I- 2.4.2 


0.753 


0.612 


Herwig-h-h 2.5.1 UE7 


0.425 


0.458 


Pythia 8 


0.777 


0.815 


Phojet 


0.643 


0.668 



and rapidity, along with the transverse momentum, lon- 
gitudinal momentum fraction, and number density as 
a function of radius for charged particles within these 
jets, using early 7 TeV LHC collision data collected with 
the ATLAS detector. The study of jets with tracks al- 
lows for precise measurements of low- momentum jets and 
their properties, thus complementing calorimeter-based 
jet measurements and allowing the study of the transition 
from soft collisions to jet production in the pcrturbative 
regime of QCD. It also provides additional observables 
for consideration in the tuning of MC event generators, 
which complement existing "minimum bias" and under- 
lying event measurements. 

No tune or model presented here agrees with all quanti- 
ties measured within their uncertainties, suggesting that 
future MC tunes may be improved. Difficulty in model- 
ing the transition between soft and perturbative physics 
is indicated by disagreements between data and all MC 
distributions in the 10-20 GeV range in the dependence 
of the charged jet cross section on jet momentum. Depen- 
dence of the cross section on rapidity is consistent with 
predictions. Particles with large transverse momentum 
pJJ?' with respect to the jet that contains them are pro- 



duced more copiously than any model predicts, as are jets 
with large charged particle multiplicity. The longitudi- 
nal momentum fraction z is best described by the Pythia 
6.421 AMBTl tune. The charged particle number den- 
sity Pch{r) is well-described by the Pythia 6.421 AMBTl, 
Pythia 8.145 4C, and Herwig-t-l- 2.5.1 UE7 tunes. With 
the exception of the Herwig-f-t- 2.4.2 default tune, which 
greatly disagrees with these measurements, all models ap- 
pear to underestimate the contribution of the underlying 
event required to model the data. 
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FIG. 3. Multiplicity of particles per charged particle jet, over the full measured rapidity range \y\ < 1.9, with anti-fct radius 
parameter R as indicated. Figure a (c) shows the distributions for five momentum ranges with R = 0.4 (0.6), and Figure b (d) 
shows the fractional difference between a range of Monte Carlo event generator predictions and the data. The distributions for 
the Perugia 2010 tune, not shown, agree qualitatively with the Perugia tune. 
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